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Two-dimensional superintegrable systems with one third order and one 
lower order integral of motion are reviewed. The fact that Hamiltonian 
systems with higher order integrals of motion are not the same in classical 
and quantum mechanics is stressed. New results on the use of classical and 
quantum third order integrals are presented in Section 5 and 6. 

1 Introduction 

The purpose of this article is to present some new results and insights on 
superintegrable systems in classical and quantum mechanics, involving 
integrals of motion that are cubic in the momenta. We also present a review 
of this field with emphasis on the differences between "quadratic" and 
"cubic" integrability and superintegrability. 
We recall that in classical mechanics a Hamiltonian system with 
Hamiltonian H and integrals of motion X a 

H = ^g ik PiPk + V(x,p), X a = f a (x,p), a = l,...,n-l (1.1) 

is called completely integrable (or Liouville integrable) if it allows n 
integrals of motion (including the Hamiltonian) that are well defined 
functions on phase space, are in involution {H,X a } p = 0, {X a ,Xb} p = 0, 
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a,b=l,...,n-l and are functionally independent ({, } p is a Poisson bracket). 
A system is superintegrable if it is integrable and allows further integrals of 
motion Y b (x,p), {H,Y b } p = 0, b=n,n+l,...,n+k that are also well defined 
functions on phase space and the integrals{i7, X±, X n -i, Y n , Y n+ k} are 
functionally independent. A system is maximally superintegrable if the set 
contains 2n-l functions, minimally superintegrable if it contains n+1 such 
integrals. The integrals Y b are not required to be in evolution with 
Xi,...X n -i, nor with each other. 

The same definitions apply in quantum mechanics but {H, X a ,Y b } are well 
defined quantum mechanical operators, assumed to form an algebraically 
independent set. 

The best known examples of (maximally) superintegrable systems are the 

Kepler-Coulomb [1,2] system V(x) = j and the harmonic oscillator 

V(x) = ar 2 [3,4]. It follows from Bertrand's theorem [5,6] that these are the 

only two spherically symmetric examples in Euclidean space. 

In both cases the integrals of motion are first or second order polynomials 

in the momenta. The Hamiltonians and integrals are the same in classical 

and quantum mechanics (after a possible symmetrization). 

A systematic search for superintegrable systems in two-dimensional 

Euclidean space E 2 was started some time ago [7,8]. Rotational symmetry 

was not imposed. The integrals were assumed to be second order 

polynomials in the momenta. Thus, the Ansatz (in quantum mechanics) was 

H=^(p) 2 + V(x) , x = (x u x 2 ) (1.2) 

2 2 

Xa=^{f?W,PiPk} + ^9iWPi + <l>aW 
i,k=l i=l 

where {, } is an anticommutator and we put 

d 

Pj = - ^^T> ^3 = x 2 pi - xip 2 (1.3) 

The commutativity requirement [H,X a ]=0 implies that X a must have the 
following form : 

X = aLj+b(L 3 p 1 +p 1 L 3 )+c(L 3 p 2 +p 2 L 3 )+d(p 2 1 -p 2 2 )+2fp 1 p 2 +(f ) (x,y), (1.4) 
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where a,...,f are constants (or X can be first order operator). 
Using transformations from the Euclidean group E(2) (they leave the form 
of the Hamiltonian (1.2) invariant) we can transform the integral (1.4) into 
one of four standard forms. If one such operator X exists then the potential 
V(xi,x 2 ) will allow the separation of variables in cartesian, polar, parabolic 
or elliptic coordinates, respectively. If two operators {X 1; X 2 }, commuting 
with H exist, we obtain four families of superintegrable potentials, [7,8], 
namely 

Vj = a(x 2 + y 2 ) + 4 + \ V n = a(x 2 + Ay 2 ) + 4 + 71/ (1-5) 
x A y z x z 

a 1 a (3 TA a l, n ,</>. . ,0,. 

V HI = - + - — —- + —- , V IV = - + -(f3cos^-)+ 1S in^-)) 

r r z 1 + cos(<p) 1 — cos(<p) r r 2 2 

Each of these potentials is multiseparable, i.e. allows the separation of 
variables in the Schrodinger equation (and in the Hamilton- Jacobi 
equation) in at least two systems of coordinates. 

A similar search for superintegrable systems in 3 dimensional Euclidean 
space E 3 gave a complete classification of all "quadratically 
superintegrable" systems in E 3 [9,10]. 

Since then many results have been obtained on superintegrable systems in 
E n for n arbitrary and for two dimensional spaces of constant curvature, so 
called Darboux spaces (2 dimensional spaces with non constant curvature 
allowing at least two killing tensors), for complex spaces, etc [11,.. .,15]. 
The overall picture that emerges is that superintegrable systems with 
quadratic (or linear) integrals of motion have the following properties that 
make them interesting from the point of view of physics and mathematics : 

1. All finite classical trajectories are closed (more generally, they are 
constrained to an n-k dimensional manifold in phase space [16]). 

2. They are all multiseparable (in classical and quantum mechanics). 

3. Their quantum energy levels are degenerate. 

4. The Schrodinger equation for quadratically superintegrable systems in E 2 
is exactly solvable and it has been conjectured that this is true in general 
[17]. The conjecture is true in all cases so far studied. 

5. The same potentials are superintegrable in classical and quantum 
mechanics. 
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6. The integrals of motion form non-Abelian algebras under commutation, 
or Poisson commutation, respectively. These algebras may be finite 
dimensional, or infinite dimensional Lie algebras. The infinite-dimensional 
ones have a special structure and it is often advantagous to view them as 
polynomial algebras, more specifically quadratic ones [18,19,20]. 
It would be interesting to establish which of the above properties hold for 
all superintegrable systems and which are restricted to Hamiltonians of the 
form (1.1) with quadratic integrals of motion. Other systems that have been 
studied include velocity dependent potentials [21-25] and particles with spin 
[26]. In this article we review the case of a scalar potential as in (1.1) but 
with third order integrals of motion. We also present some new results. 

2 Integrability with third order integral of 
motion. 

In 1935 J. Drach published two articles on two-dimensional Hamiltonian 
systems with third order integrals of motion [27]. The main feature of these 
papers are 

1. The author found 10 such integrable potentials, each depending on 3 
constants, but not involving any arbitrary functions. 

2. The results were obtained in classical mechanics and are in general not 
valid in quantum mechanics. 

3. The space considered was a complex Euclidean space E 2 (C). 

4. The articles are very short and "condensed". Basically they are research 
announcements, but no detailed follow up articles were published. 

5. Some assumptions were made in the calculation so it is not clear whether 
the obtained list is complete. 

Much more recently it was shown by Rahada [28] and Tsiganov [29] that 7 
of the 10 Drach potentials are actually "reducible". The 7 corresponding 
Hamiltonian systems are superintegrable, in that they allow two second 
order integrals. Their Poisson commutator is the third order integral found 
by Drach. 

Since then many examples of systems with third and higher order integrals 
have been published, mainly in classical mechanics[30,...,34]. 
In 1998 J.Hietarinta published a remarkable article "Pure quantum 
integrability" [35] in which he showed that potentials exist that are 
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integrable in quantum mechanics but have free motion as their classical 
limit. The integrals of motion are third or higher order polynomials in the 
momenta and the potentials have the form V(x, y) = h 2 f(x, y) where K is 
the Planck constant and f(x,y) remains finite in the limit H — > 0. 
A systematic search for integrable and superintegrable systems in both 
classical and quantum mechanics was initiated in Ref 36. The problem was 
formulated in quantum mechanics, the classical case was considered as a 
limit for H — > 0. The Hamiltonian and first integral have the form 

H=^(pl+pl) + V(x 1 ,x 2 ) (2.1) 



X = -J2fj^i,^)p{p k 2 0<j + k<3 (2.2) 

where pj are as in (1.3) and the potential V and functions fjk are to be 
determined. The commutativity requirement [H,X]=0 implies that for any 
potential the third order integral X will have the form 

X= A ijk{Ll,p[p2} + {9i(xi,x 2 ),Pi} + {g 2 (x 1 ,x 2 ),p 2 } (2.3) 

i+jr' + fc=3 

where {, } is an anticommutator, L 3 is the angular momentum, as in eq. 
(1.3). The constants and functions V, <?i and g 2 are subject to further 
determining equations : 

(gi)x 1 = 3fi(x 2 )V Xl + f 2 (xi,x 2 )V X2 , (g 2 ) X2 = f3(xi,x 2 )V Xl + 3/4(^1)^ 

(2.4) 

(9i)x2 + (92) X1 = 2(f2(xi,x 2 )V xl + f 3 (x 1 ,x 2 )V X2 

and 

h 2 

9iV Xl + g 2 V X2 = —(fiV XlXlXl 

foMx\X\X2 f 3^X1X2X2 f 4^X2X2X2 ) (2-5) 

+8A 300 ( Xl V X2 - x 2 V Xl ) + 2(^210^ + A 201 V X2 ) 
The functions are defined in terms of the constants A^ k as 

fi = -^300^2 + ^210^2 - A 120 x 2 + A 030 (2.6) 
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h = 

h = - 



■■ 3A 300 ZiX2 - 2A 2W X 1 X 2 + A 201 xj + ^120^1 ~ ^111^2 + A)21 

-3^300X^X2 + A 2W x\ - 2A 2Ql xix 2 + Amxx - A i02 x 2 + A 0l2 

Ia = ^300^1 + ^210^1 + A W2 Xi + A 003 



We can sum up the results on third order integrability as follows. 

1. The leading part of the integral X lies in the enveloping algebra of the 
Euclidean Lie algebra e(2). 

2. Even and odd terms in X commute with H separately. We required the 
existence of a third order integral. A second order one of the form (1.4) 
may, or may not exist. 

3. The 10 constants and 3 functions gi,g 2 and V are to be determined from 
the (overdetermined) systems of 4 equations (2.4) and (2.5). The fact that 
the system is overdetermined implies the V(x,y) will satisfy some 
compatibility conditions. Indeed, one can obtain 1 linear third order 
equation for V and 3 nonlinear ones [36]. The existence of a third order 
integral is hence much more constraining than the existence of a second 
order one ( in agreement with Drach's results). 

4. The classical and quantum cases differ ! Indeed, the determining equation 
(2.4) are the same in both cases, but eq. (2.5) involves the Planck constant 
h. In the classical limit h — > it simplifies greatly ( the right hand side goes 
to zero). 

5. Instead of attempting to find a general solution we turn to a simpler 
problem, namely that of superintegrability. We shall consider potentials 
V(x,y) for which two integrals of motion exist, one third order one and one 
that is either a first, or second order polynomial in the momenta. 



3 Systems with one first order and one third 
order integral 

Let us now assume that the Hamiltonian (2.1) allows a first order integral 
of motion Y. This means that the potential has a geometric symmetry : it 
is invariant either under rotations, or translations in one direction. With no 
loss of generality we can consider just two representative cases : 



6 



1) V = V(r), Y = L 3 , r=^x 2 + y 2 , (3.1) 

2) V = V(x), Y = P 2 . 



Eq. (2.4) and (2.5) for the existence of a third order integral X greatly 
simplify in both cases. 

Let us first consider classical mechanics, i.e. H = in eq. (2.5). Solving the 
systems (2.4) and (2.5) for the above potentials, we obtain 

f\l ft 

V 1 = -, V 2 = ar 2 , V 3 = ax, V 4 = — . (3.2) 

All of these potential are well know [1,...,8,37] as quadratically 
superintegrable. Indeed, the third order integral in all these cases is the 
Poisson commutator of a second order integral with the square of the first 
order one. 

In quantum mechanics we again obtain the potentials (3.2) but in addition 
we obtain a potential V(x) satisfying 

h 2 {V"{x)) 2 = AV 3 + aV 2 + f3V + 1 = 4(V - A 1 )(V - A 2 )(V - A 3 ) (3.3) 

where a,j3 and 7 (and Ai,A 2 ,A 3 ) are constants. Eq. (3.3) is solved in terms 
of elliptic functions. Depending on properties of the roots Ai and on the 
initial conditions, we obtain 3 types of solutions : 

K=(IMV m W), V- 2(m m ) + l) (3-4) 

v = (hu) 2 ! 



sn 2 (ux, k) 

where sn(ux, k) and cn(ux, k) are Jacobi elliptic functions. If two of the 
roots Ai coincide, we get potentials expressed in terms of elementary 
functions 

y_ (M 2 y _ (^) 2 (3 5) 

(cosh{uox)) 2 ' (sinh(ux)) 2 
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We see explicitly that for H — > all these potentials vanish. The potentials 
(3.4) and (3.5) are "irreducible" and genuinely superintegrable. The two 
algebraically independent integrals are 

y = P 2 , X = {L 3 ,p 2 1 } + {(a-3V)y, Pl } + {-ax + 2xV + J V(x)dx,p 2 }, 

(3.6) 

a = A 1 + A 2 + A 3 

Even though the potentials V(x) in (3.4) and (3.5) are one-dimensional 
(depend on x alone), their superintegrability is a two-dimensional 
phenomenon, since the integral X involves rotations L%. 



4 Systems with one second and one third 
order integral 

Let us now assume that the Hamiltonian (2.1) allows one second order 
integral of motion Y. This integral will have the form (1.4). As mentioned 
in the Introduction, 4 classes of such Hamiltonians exist [7,8]. Here we shall 
restrict to one of the 4 classes, namely to potentials allowing the separation 
of variables in cartesian coordinates. Thus we have 

y=\{p\-pl) + V 1 {x)-V 2 {y), V(x,y) = V 1 (x) + V 2 (y). (4.1) 

We substitute the above potential into the determining equations (2.4) and 
(2.5) in order to determine potentials that also allow a third order integral. 
The obtained system of equations is quite manageable and was analyzed by 
S. Gravel [38]. 

The results in the classical and quantum cases are very different. To 
illustrate the difference, let first consider an example. One particular 
solution of the determining system (2.4), (2.5) is the potential 

V(x,y) = ^-y 2 + V(x) (4.2) 
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where V(x) satisfies a fourth order nonlinear ordinary differential equation 

H 2 V"" = 12u 2 xV + Q(V 2 )" - 2u 2 x 2 V" + 2oo 2 x 2 . (4.3) 

We mention that this equation is also obtained as a nonclassical reduction 
of the Bousinesq equation [39, 40]. For h ^ this equation is solved in terms 
of the Painleve transcendent [41] Piv{x, — ^j-)- The classical limit is 
singular (H — > reduces the order of the equation ). In the classical case we 
can reduce to quadratures and obtain a quartic algebraic equation for V(x) : 

- W\x) + Uu 2 x 2 V 3 (x) + (Qd - ^-cu 4 x*)V 2 (x)+ (4.4) 

6^,6 _ 2duJ 2 x 2 )v{x) + cx 2 _ d 2_ ^4 _ 1 8^,8 = q 

2 2 16 

For special values of the constants a,...,d,u; namely 

8oo s b 3 , uj A b 2 , , 

c=^~, d=^- (4.5) 

eq. (4.4) has a double root and we obtain 

2 2 2u 

V li2 = ^(2b + 5x 2 ±AxVbTx^), V 3 = V, = ^-x 2 -^ r (4.6) 

A complete analysis of eq. (2.4) and (2.5) is given in Ref. 38. All together 8 
classical superintegrable systems that are separable in cartesian coordinates 
and allow at least one third order integral of motion exist. Three of them 
are reducible (the harmonic oscillator and potentials Vj and Vn of eq. (1.5)) 

The 5 irreducible classical potentials are 

y=^(9x 2 + i/ 2 ) 

V = [3 2 yf\x~\+f3 2 v^\ 

V = a 2 \y\ + fr 2 v4x| 

V = ^y 2 + V(x) 

V = a\y\ + f(x), 

where V(x) satisfies eq. (4.4) and f(x) satisfies f 3 — 2bxf 2 + b 2 x 4 f — d = 0. 
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The quantum case is much richer : 21 superintegrable cases of the 
considered type exist, 13 of them irreducible. The potentials are expressed 
in terms of rational functions in 6 cases, elliptic functions in 2 cases and 
Painleve transcendents [41] Pi,Pn and Piy in 5 cases. Let us just present 
the irreducible potentials. 

Rational function potentials : 

v — h 2 \ x2+ y 2 4- 1 i i i 

V = r?[&(a? + y 2 ) + £ + j^p + (^] 

V = h i~8^( X + V ) + fe+aF + (F^F + {x+a) 2 + {x-a) 2 \ 

V=^(9x 2 + y 2 ) 

V = ^(9x 2 + y 2 ) + f 

y _ fc2rtoW _j l , 3_l 

V - n L 8a 4 + (y-a) 2 + (y+a)2 J 

Elliptic function potentials : 

V = h 2 P(y) + V(x), V(x) is arbitrary 

F = h 2 (P(y) + P(x)), P(x) is a Weierstrass elliptic function 

Painleve transcendent potentials : 

V = ^(ufPj^x) + U 2 Pj(uJ2y)) 

V = ay + ffufP^Uix) 

V = bx + ay+(2hb)IP 2 I ((§)^x) 

V = ay+ {2h 2 b 2 )HP' H {{^)hx) + 

V = a(x 2 + y 2 ) + f {P' IV {^x) - l^Pf y (^x) - l^ x p IV (=^ x )) 

The occurence of Painleve transcendents as superintegrable potentials 
seems somewhat surprising. It is less so once we remember the relation 
between the Schrodinger equation and the Korteweg-de Vries equation [42]. 
Solutions of the KdV include Painleve transcendents. Painleve transcendent 
potentials have already occured in other contexts [43,44,45]. 

5 Trajectories for the classical systems 

The trajectories for all 8 superintegrable classical potentials were presented 
in Ref. 46. The trajectories can be obtained by integrating the equations of 
motion and then eliminating the time variable from the result. 
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Alternatively, the trajectories can be obtained directly from the integrals of 
motion. Indeed, in all cases we have 



-p\ + f{x) = E 1) -pl + g( x ) = E 2 



(5.1) 



X = fJ.pl + vp\v2 + PPiPl + crpl + #1 + #2 = K 



(5.2) 



where //,... ,<7 are known low order polynomials in x and y, f(x),g(x),0(x, y) 
and ip(x,y) are known functions and E±,E 2 and K are constants depending 
on the initial conditions. 

From eq. (5.1) we express P\ and P 2 in terms of Ei,E 2 ,i(x) and g(y) and 
substitute into eq (5.2). This directly gives us an equation for the trajectory 
(not necessarily in a convenient form). If the functions f(x) and g(y) satisfy 
f(x) > 0, g(y) > 0, for x 2 > Xo, y 2 > yo for some constants Xq and yo, then 
the motion will be bounded. 

From the figures of Ref 46 we see that the finite trajectories are all closed, 
as predicted by the general theory. 

6 Algebras of the quantum integrals of motion 

The integrals of motion in quantum mechanics form associative polynomial 
(cubic) algebras with respect to Lie commutation. The commutator 
relations in all cases are of the form 



[A,B]=C, [A,C]=aB, [B, C] = {JA 3 + 1 A 2 + 5A + e (6.1) 



where 7, 6 and e are polynomials in H (for the classical cubic algebras of 
Poisson brackets see Ref. 46). 

Inspired by the work of Daskaloyannis and al. on quadratic algebras 
[18,47,48] we construct a realization of the cubic algebra (6.1) using a 
deformed oscillator algebra and Fock space basis {6*, 6, N}. We put 



[N, 6*] = b\ [N, b] = -b, b l b = $(iV), bb l = <S>{N + 1) (6.2) 



where $(iV) is the so called "structure function". To construct 
representations of the cubic algebra we impose conditions on $(7V). We 
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require that a natural number p should exist, such that 

$(p+l) = (6.3) 

and further we impose 

$(0) = 0, > 0, x > (6.4) 

We put 

A — A(N), B = b(N) + b t p(N) + p(N)b (6.5) 

where the functions A(N),b(N),p(iV) and <3>(iV) are to be determined from 
the cubic algebra (6.1). 

The next step is to construct a finite dimensional representation of this 
algebra with A, N and K (the Casimir operator) diagonal. We obtain a 
realization of a parafermionic oscillator 

N\k, n >= n\k, n >, K\k,n >= k\k,n > 

A\k,n >= A(k,n)\k,n > 

A\k,n>= VS(n + u)\k,n> (6.6) 

$(0,u,Jfc) = 0, $(p+ l,u, k) = 

A complete discussion of this approach to quantum superintegrable systems 

with cubic integrals of motion will be presented elsewhere. 

Here we just present one example showing how one can obtain energy 

spectra. Let us consider one of the rational quantum potentials of Section 4, 

namely 

v = h 2 { ^W _1 ] (67) 

8a 4 (x-a) 2 (x + a) 2 v ' 

The Casimir operator K and structure function <&(#) are in this case found 
to be 

h 4 h 6 h 8 h 10 

K = -16h 2 H 4 + 32^H 5 + 16^jH 2 -40^-H-3^ . (6.8) 

a 2 a 4 a 6 a 8 
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(6.9) 

The condition $(0, u, k) — is now an equation for u and $(p + 1, u, k) — 
provides the energy spectrum. In the case of the potential (6.7) the 
constant a can be real, or pure imaginary (V(x,y) is real in both cases). The 
energy spectrum we obtain is 

ia , a G R, pGN 
a G R, peN 

7 Conclusion 

The main conclusion that we draw at this stage is that going beyond 
quadratic integrability for scalar potentials opens new horizons and poses 
new questions. 

One of the questions is : What does one do with integrals of motion that do 
not lead to the separation of variable in the equations of motion ? A partial 
answer is given in this article and is very different in classical and quantum 
mechanics. We have seen that the integrals can be used to calculate 
trajectories directly without solving the classical equation of motion. In 
quantum mechanics one uses the algebra of integrals of motion to obtain 
information about energy spectra and wave functions. 

Acknowledgments The research of P.W. is partially supported by a 
research grant from NSERC of Canada. This article was written while he 
was on sabbatical leave at ENEA, Frascati, Italy and Universita di Roma 
Tre. He thanks G.Dattoli and D.Levi for hospitality. The research of I.M. is 
supported by a doctoral research scholarship from FQRNT of Quebec. This 
article was written while he was visiting the Universita di Roma Tre. He 
thanks D.Levi for his hospitality. 



E = 



Hp + 2) 



2a 2 



a = 



E 



Hp + 3) 
2a 2 



13 



8 References 



1. V.Fock, Z.Phys. 98, 145-154 (1935). 

2. V.Bargmann, Z.Phys. 99, 576-582 (1936). 

3. J.M.Jauch and E.L.Hill, Phys.Rev. 57, 641-645 (1940). 

4. M.Moshinsky and Yu.F.Smirnov, The Harmonic Oscillator In Modern 
Physics, (Harwood, Amsterdam, 1966). 

5. J.Bertrand, C.R. Acad. Sci. Ill, 77, 849-853 (1873). 

6. H.Goldstein, Classical Mechanics, (Addison- Wesley, Reading, MA, 1990) 

7. J.Fris, V.Mandrosov, Ya.A.Smorodinsky, M.Uhlir and P.Winternitz, 
Phys.Lett. 16, 354-356 (1965). 

8. P.Winternitz, Ya.A.Smorodinsky, M.Uhlir and I.Fris, Yad.Fiz. 4, 625-635 
(1966). (English translation in Sov. J.Nucl.Phys. 4, 444-450 (1967)). 

9. A.Makarov, Kh. Valiev, Ya.A.Smorodinsky and P.Winternitz, Nuovo 
Cim. A52, 1061-1084 (1967). 

10. N.W.Evans, Phys.Rev. A41, 5666-5676 (1990), J.Math.Phys. 32, 
3369-3375 (1991). 

11. E.G.Kalnins, J.M.Kress, W.Miller Jr and P.Winternitz, J.Math.Phys. 
44(12) 5811-5848 (2003). 

12. E.G.Kalnins, W.Miller Jr and G.S.Pogosyan, J.Math.Phys. A34, 
4705-4720 (2001). 

13. E.G.Kalnins, J.M.Kress and W.Miller Jr, J.Math.Phys. 46, 053509 

(2005) , 46, 053510 (2005), 46, 103507 (2005), 47, 043514 (2006), 47, 043514 

(2006) . 

14. E.G.Kalnins, W.Miller Jr and G.S.Pogosyan, J.Math.Phys. 47, 
033502.1-30 (2006), 48, 023503.1-20 (2007). 

15. M.A.Rodriguez and P.Winternitz, J.Math.Phys. 43(3), 1309-1322 (2002). 

16. N.N.Nekhoroshev, Trans. Moscow Math.Soc 26, 180 (1972). 

17. P.Tempesta, A.V.Turbiner and P.Winternitz, J.Math.Phys. 42, 4248 
(2001). 

18. C.Daskaloyannis, J.Math.Phys. 42, 110 (2001). 

19. Ya.I.Granovskii, A.S.Zhedanov and I.M.Lutzenko, J.Phys. A24, 3887 
(1991). 

20. P.Letourneau and L.Vinet, Ann.Phys. 243, 144 (1995). 

21. B.Dorizzi, B.Grammaticos, A.Ramani and P.Winternitz, J.Math.Phys. 
26, 3070-3079 (1985). 

22. E.McSween and P.Winternitz, J.Math.Phys. 41, 2957-2967 (2000). 
23 J.Berube and P.Winternitz, J.Math.Phys. 45, 1959-1973 (2004). 



14 



24. F.Charest, C.Hudon and P.Winternitz, J.Math.Phys. 48, 012105.1-16 
(2007). 

25. G.Pucacco and K.Rosquist, J.Math.Phys. 46, 012701 (2005). 

26. P.Winternitz and I.Yurdusen, J.Math.Phys. 47(10), 103509 (2006) and 
ArXiv math-ph. :0711.0753 (2007). 

27. J.Drach, C.R.Acad.Sci.III, 200, 22 (1935), 200, 599 (1935). 

28. M.F.Ranada, J.Math.Phys. 38, 4165-4178 (1997). 

29. A.V.Tsiganov, J.Math.Phys. A33, 7407 (2000). 

30. G.R.Holt, J.Math.Phys. 23(6), 1037-1046 (1982). 

31. G.Pucacco and Rosquist, J.Math.Phys. 46, 052902.1-21 (2005). 

32. A.S.Fokas and P.A.Langerstrom, J.Math.Anal.Appl. 74, 325-341 (1980), 
74, 342-358 (1980). 

33. J.Hietarinta, Physics Reports 147, 87-154 (1987). 

34. G.Thompson, J.Math.Phys. 25(12), 3474-3478 (1984). 

35. J.Hietarinta, Phys.Lett. A246, 97-104 (1998), J.Math.Phys. 25, 
1833-1840 (1984). 

36. S. Gravel and P.Winternitz, J.Math.Phys. 43(12), 5902 (2002). 

37. M.B. Sheftel, PTempesta and P.Winternitz, J.Math.Phys. 42(2), 
659-673 (2001). 

38. S. Gravel, J.Math.Phys. 45(3), 1003-1019 (2004). 

39. P.A.Clarkson and M.D.Kruskal, J.Math.Phys. 38(11), 5944-5958 (2007). 

40. D.Levi and P.Winternitz, J.Phys. A22, 2915-2924 (1989). 

41. E.L.Ince, Ordinary Differential Equations (Dover, New York, 1944). 

42. M.J.Ablowitz and P.A.Clarkson, Solitons, Nonlinear Evolution 
Equations and Inverse Scattering, (Cambridge Univ. Press, 1991). 

43. W.I.Fushchych and A.G.Nikitin, J.Math.Phys. 38(11) 5944-5959 (1997). 

44. A.P.Veselov and A.Shabat,Funkt.Analiz.Prilozh 27(2), 1-21 (1993). 

45. A.P.Veselov, J.Phys. A34, 3511 (2001). 

46. I.Marquette and P.Winternitz, J.Math.Phys. 48(1) 012902 (2007). 

47. D.Bonatsos, C.Daskaloyannis and K.Kokkotas, Phys.Rev. A48, 
R3407-R3410 (1993). 

48. C.Daskaloyannis and K.Ypsilantis, J.Math.Phys 47, 042904.1-38 (2006). 



15 



